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A CLASSICAL  MODEL  OF  FERROMAGNETISM  AND  ITS  SUBSEQUENT  QUANTIZING 
IN  THE  AREA  OF  LOW  TEMPERATURES 

Article  by  G.  Heller  and  H.  A.  Kramers,  Utrecht 
— Communicated  at  the  Meeting  of  June  30,  1934. 

AMSTERDAM  ROYAL  ACAD.  SCI.,  Vol  37:378-385(1934). 

1.  Summary:  The  classical  model  of  ferromagnetism  is  defined. 
Only  in  the  case  of  the  space  lattice,  it  clearly  results  in 
ferromagnetic  saturation.  The  simple  quantizing  of  this  model 
results  exactly  in  the  formulas  of  Bloch*  and  Miller**. 

* Bloch,  ZS.  F.  PHYSIK,  Vol  61:206(1930)  and"  Vol  74:295(1932). 
**  Miller,  ZS.  F.  PHYSIK,  Vol  82:559(1933). 

2.  Classical  Theory  of  Ferromagnetism,  Especially  at  Low 
Temperatures 

We  consider  a lattice  of  N atoms.  Let  be  the  aperator  of  the 

mechanical  moment,  / e h V 

— inc  2 ) the  operator  of  the  magnetic 

moment  of  the  i-th  atom.  If  H is  the  external  magnetic  field 

strength,  then  the  energy  operator,  as  far  as  it  is  function  of 

the  spin  moments: 

’ w-  — uH  £ o,  — i J /*(®»  «»)••••  0) 

isi  »<• 


The  I are  exchange  integrals  between  the  atoms  s and  t, 


If  the  eigenvalues  E of  this  operator  are  known,  the  condition 

z 

sum  - 

Z[6.H)=  £e~~i  [e  = kT)  . .....  (2) 

can  be  calculated  and  the  expected  value  of  the  magnetic  moment 


* r*6  ■■ 


is  found. 


With  Heisenberg  and  Bloch,  we  assume  that  the  I are  different 

s t 

only  for  adjacent  units  of  zero  and  we  establish  the  same  task  as 
Bloch,  — the  calculation  of  p for  the  linear  chain,  the  square 
layer  lattice  and  the  cubis  space  lattice. 

The  "classical"  treatment  of  the  calculation  consists  of  the  fact 

that  we  consider  the  operators  of  equation  (1)  as  c-numbers  in 

the  sense  of  classical  mechanics.  The  vector  s.(|s|=*n)  with  the 

2 2 2 1 

directional  cosines  X^Y^Z^(X  +Y  +Z  =1)  corresponds  to  the  spin 
operator  o4  and  equation  (1)  is  for  the  linear  chain: 

W-  — /tnH  2 Zk  — iln*  2 (*»  *»♦,  + Yk  YM  + Z*  Z*,)  (4) 

knl  *=l 

if  we  select  the  field  direction  as  Z direction,  whereby  we  set 
Xjj+^“X^  and  so  forth  for  the  completion  of  the  lattice. 

_W 

We  now  obtain  the  condition  3uro  through  the  integration  of  e-0 

through  the  canonical  coordinates  of  the  entire  system.  If,  for 

example,  we  describe  the  spin  orientation  of  the  i-th  atom  by  the 

canonical  coordinates  A 

2*a  and^?^  and  introduce  the  differential 

of  the  space  angle  by 

n dm  = dsn  dm 

we  obtain  the  condition  sum: 

Z(0.//)  = ( J « * dwt  Jci>, du>N.  . . . (5)| 

The  exact  evaluation  of  (5)  presents  considerable  mathematical 
difficulties.  We  limit  ourselves  to  the  investigation  of  the 
behavior  at  low  temperatures  6<<I. 
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We  can  then  assume  that  almost  all  spins  lie  almost  in  the  Z 
direction.  We  can  then  develop  the  energy  W by  powers  of 
Y ^ and  we  obtain  with 


and 


W=-Nn(jj  + f‘H)  + I-Y  [Zy+°)(Xl+Yh-X>X>*x-Y>  MW 
dui,  — dXi  dYi 

jgr  — J' e~TdXt...dXNdYt...dYs‘ 


In  practice,  the  integrand  is  markedly  different  from  0 only  when 

% 

all  X^,  have  very  small  amounts.  We  may  therefore  extend  the 
integration  limits  from  to  +®  (This  conclusion  is  of  course 
only  permissible  when  the  integral  remains  convergent). 

We  apply  the  integral  formula: 

J t~Xm+xt  *»£/*,  dX7...  dX»  = ||A||-V.s=^  j/^k  . (7) 

In  it,  are  the  coefficients  of  a symmetrical  matrix,  | A | is 

its  determinant,  A^(A=1...N)  are  its  eigenvalues. 

Proof:  There  are  real  unitary  matrices  U,  which  transform  A on 

principal  axes: 

U~l  AU—A.  \ 

If  we  introduce  new  coordinates  through: 

x>  = 2 U»  Xi 


we  obtain: 
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Hith  (7),  the  condition  sum  becomes: 


hd  Y_j£L 

Z = e»^  -jO~ 


W! 


The  A^  are  the  eigenvalues  of  the  square  form  within  brackets  in 
(6).  They  are  found  from  the  secular  equations: 

— iu»-u  + (l  + a—jl;)im“i  =®  • • • • (9) 

i 

One  solution  is: 


*-i, : 


Uti  = eN'  ,At  — l-\-a — cos  jj  l 


(10) 


Equation  (10)  provides  us  with  N independent  eigenvectors,  thus 

a complete  system.  We  thus  also  know  all  eigenvalues  and  thus, 

* 

according  to  (8),  the  condition  sum: 

(V tt  “ i h (' + «— ■**)• 

In  the  case  of  large  numbers  of  atoms  N,  we  can  replace  the  sum 
by  an  integral  through 

- *» 

h Z=  ^ (y  + ,.//) + Nl«  ^ - gj'la (1  -f-.-co. £) di  . (II), 

0 

According  to  (3)  and  (11) , the  expected  value  of  the  magnetic 
moment  becomes: 


- dtnZ  .»  N f*6  f d± 

= n/J  l+o- 


cos  ( 


= Nnn-N^ 


1 


«/  |/o  (a + 2) 


This  expression  diverges  with  respect  to  -•  for  H-*-0,  i.e.  a+0 . 

This  result  does  of  course  not  have  any  physical  significance,  but 
it  perhaps  makes  possible  the  conclusion  that  there  is  no  tenden- 
cy to  ferromagnetism. 


■L 
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Analogous  calculations  apply  to  square  layer  and  space  lattices. 
The  energy  of  the  simple  square  layer  lattice  is: 


W- 


=■ — n/*H  Z Zki  — -Ln*  Z (XuXma  + XkjX».Hi~h 


2 I.M 

+ Yu  YM.i+  Yu  Yk.1+1  + Zu  Zku.  / +Zu  Zt,.  i-bt) 


With  the  same  neglected  terms  as  above,  the  condition  sum  becomes: 

«N’ 


Z=e 


\J  n) 


The  secular  problem  for  the  determination  of  the  A^y  is: 

% 

i(llk-l.t4«  + (2  + « — Al^Uk.l,,*  —\(tlk*U;>r  + = 

The  solutions  are: 


BU»=f 


■dip  — 2 + a — cos  ~ kl—cos  ~ I ft 


N 


We  set: 


2.t  t 2* . 

A ^ */*■“  *1 


and  obtain  for  the  logarithm  of  the  condition  sum: 
lnZ  — ^~{nl+ftH)  + N2ln~Y  — 


-ms 


ln(2  + a — cosk  — cos*i)d$dt) 


The  expected  value  for  the  magnetic  moment  thus  becomes  for  the 
layer  lattice: 


/*  = N'm*  ~ (SJ^T //  2 + o -~cos(— i 


and  in  an  analogous  manner  for  the  space  lattice: 
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- xn  (N\'nB  CCC  dtdtidC 

P — Wnn—  „I  JJJ  3+a-co»$  — cost -cost 

For  a-»-0,  only  the  space  lattice  retains  a finite  saturation: 


— (1  — a0)  i 


. 02) 


while  the  quantum-mechanical  calculation  resulted  in 


— — a*  pi.) 


This  calculation  shows  that  the  ferromagnetism  is  not  a typical 
quantum-mechanical  effect.  While  the  magnetic  saturation  at  low 
temperatures  is  considerably  enhanced  by  quantum  effects,  it  also 
occurs  in  the  "classical"  consideration,  under  the  same  condition 
as  in  the  Bloch  calculation. 


As  already  noted  by  Bloch,  (12)  contradicts  the  Nernst  law  of 
thermodynamics.  It  is  reasonable  to  provide  relief  through  subse 
quent  quantizing,  in  a process  analogous  to  the  Debye  theory  of 
specific  heat. 

3.  The  Quantizing  of  Spin  Waves 


With  small  deviations  and  m. <>«... 

field  direction,  the  coordinates  = 
conjugated  canonically.  We  obtain  a new  system  of  canonical  va- 
riables in  that  we  subject  p^  and  q^  to  contragredient  transforma 
tions.  In  accordance  with  (7),  we  select  U real  and  form 

Pk  — 2 “**  P\  9k  — 2 (U-'U  q\  = 2 an  q\ 


of  the  maenetic  moments  from  the 

are 


and  q^U^-Y, 


In  the  eliminated  variables,  we  obtain  for  the  energy  of  the 
linear  chain: 

W=  ~ Nn  (y  + *")  + J ,|  (»  +«  -cos^l)  {pf  + q?) 

As  this  equation  shows,  our  system  is  a harmonic  oscillator  in 
each  of  the  N pairs  of  variables  p'A,q'x,  namely  with  the 
frequency : 

2/n  /,  , 2n  A 

an  = ^-[i+a-cos—iJ (13) 

Each  magnetic  condition  of  the  lattice  near  the  absolute  satura- 
tion can  be  considered  as  superposition  of  "spin  waves"  of  these 
fundamental  frequencies. 

An  eigenvalue  of  the  entire  system  is  obtained*  by  ascribing  an 

arbitrary  quantum  number  n to  each  oscillator  u>  : 

A X 

— Nn^‘2+<*Hj  + lnihci>i.  . . . (14) 

Equation  (14)  contains  all  eigenvalues  of  the  system.  According 
we  calculate  the  condition  sum: 

- - - _ 

«.=•  «N=o 

According  to  (14),  we  have: 

Z-t  • 17  2e  * 

i*=l  Hj=4 

»•&*'")  N , 

= e * II U— 

i*i  hui 

1— e~  * 

hZ=tr(j+l‘Hyi'"  "'1) 
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We  introduce  the  according  the  (13)  and  again  replace  the  sum 
by  an  integral: 


The  expected  value  of  the  magnetic  moment  for  the  linear  chain 
now  becomes: 


(£)**£ — 


The  analogous  calculation  results  in: 


for  the  layer  lattice  and 


cr 

d£d>) 

JJ 

0 « 

j 

I.T 

rr 

dSdtjdi 

0 • 


for  the  space  lattice. 


This  is  exactly  the  formula  of  Miller  if  the  integral  is  still 
evaluated  by  the  series 


cos  { = 


i 


In  this  connection,  it  is  perhaps  not  uninteresting  to  note  that 
the  Miller  calculation  can  be  interpreted  in  such  a way  that  only 
those  condition  of  the  lattice  are  taken  into  consideration  in 
which  almost  all  atoms  in  the  direction  of  the  magnetic  field 
possess  the  spin  moment  JL~h 

(v  2 and  the  few  remaining  ones  the  spin 

— — 1 \ ft ' 

2 J , thus  only  less  by  one  quantum  number.  In  the 
classical  picture,  the  assumption  thus  corresponds  that  the  spin 
directions  deviate  only  little  from  the  field  direction. 


-9- 


The  complete  agreement  of  the  above  calculation  with  the  theory 

of  Bloch  results  in  a very  clear  interpretation  of  the  Bloch  spin 

waves.  Their  carrier  is  a system  of  spins  elastically  coupled  by 

exchange  forces.  An  excitation  of  the  n-th  quantum  condition  of 

the  classically  calculated  eigenf requency 

2 In  f 2 » 2 n 2a  \ : 

Wiri#.  — •*=“  f 3 + a — COJ * — COS  * COS  jy/tJ.j 

corresponds  to  the  condition  that  n spin  deviations  according  to 
the  De  Broglie  wave 

eN 

wander  through  the  lattice. 

* 

I 

The  Bose  statistic  of  the  elastic  natural  oscillations  is  identi- 
cal with  the  Bloch  "sum  of  Bose  gases". 

\ 

However,  it  should  be  noted  that  the  exact  agreement  of  our  cal- 
culation with  the  results  of  Bloch  and  Miller  was  hardly  clear 
from  the  start  because  the  normal  coordinates  p'^,q'  do  not  at 
all  correspond  to  well-defined  quantum-mechanical  operators. 

This  is  probably  due  to  the  fact  that  we  had  to  omit  the  null 
point  energy  in  (14)  in  order  to  arrive  at  agreement  with  the 
quantum-mechanical  calculation  with  regard  to  the  absolute  value 
of  the  energy  and  of  the  magnetic  saturation. 


